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Abstract 

Short-distance expansion of the total semileptonic B widths is reviewed for the OPE- 
conformable scheme employing low-scale running quark masses. The third- and fourth- 
order BLM corrections are given and the complete resummation of the BLM series pre- 
sented. The effect of higher perturbative orders with running quark masses is found very 
small. Numerical consequences for \Vcb\ are addressed. 
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Total decay rates of heavy flavors are dominated by short-distance physics; nonpertur- 
bative effects originating from a typical hadronic momentum scale are described by the 
OPE. The semileptonic B widths Tsi{b^c) and Tgi^b^u) offer an unmatched accuracy in 
determination of the corresponding CKM elements. To make a full use of the infrared-free 
nature of perturbative corrections to the widths they should be computed in the context 
of the Wilsonian OPE; in particular this suggests using the running quark masses nor- 
malized at a low scale fi<^mii [1]. The perturbative corrections (i.e., the ag-expansion 
of the Wilson coefficients of the unit operator QQ) explicitly depend in this approach on 
the separation (normalization) scale /i. 

The //-dependence of the Wilson coefficients is uniquely determined [2, 3] by normal- 
ization point independence of the inclusive width and by known //-dependence of quark 
masses and local heavy quark operators in the OPE series. Below the expressions are col- 
lected for the scheme based on running 'kinetic' heavy quark masses and other operators 
defined via the small velocity (S V) sum rules [4, 5] . These definitions are complete even 
beyond perturbation theory; likewise they are valid to any perturbative order. 

1 Perturbative corrections 

Considering pure short-distance (perturbative) effects, we denote 

^Tib^c) = \Vcb\' (1 + Aw) {zo{r) + a,(r;/x)^ + a2(r;/.) 

+ a3(r;/.)(^)%...} , (1) 

where ZQ{r) is the tree-level phase space factor and r — mKii) /mKii): 

Zo{r) = 1 - 8r + 8r^ - - 12r^ In r . (2) 

In what follows is the MS strong coupling normalized at m^, unless indicated explicitly. 

The electroweak corrections A^^ describe the ultraviolet renormalization of the semilep- 
tonic Fermi interaction and, being the QED-counterpart of the familiar factors c± in 
nonleptonic Lagrangian, has a well-known logarithm of Mz/mh due to virtual photon 
exchanges. It depends only on the fermion charges and is the same as in /9-decays of 
neutron or hyperons [6]: 
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The /f-dependence of the perturbative factors ai(r; /i) and a2(r; /i) through the second 
order in ag is given by 

ai(r;/x) = af\r) + [5Ai + 3p^] Zo{r) + [2(VF-r)Ai + {l-r)p,] (4) 



a2{r;ii) = 4'^(r) + [5Ai + 3pi] ar{r) + 2(V^-r)Ai + (l-r)pi 
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where a^k\r) = afc(r;0) are the coefficients perturbatively computed at n = 0, i.e. for the 
pole quark masses and the pole-type definition of the kinetic and higher-order operators. 
The b^u case corresponding to r = amounts to using zo — 1, ai — — | (tt^ — ^) while 
discarding all terms with derivatives. 

Four parameters Ai^2 and pi^2 entering above denote the coefficients in the perturbative 
expansion of A(/i)/mb and ii'^{iJ,)/ml to first and second order in as{M)/Tr, respectively. 
In the the kinetic mass scheme they are [7] 
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Here 



The general structure of Eqs. (4) and (5) holds in higher orders in ag, but the exphcit 
coefficients change and the terms proliferate. This does not happen for the BLM cor- 
rections^ to any order, for which the translation between different fi takes basically the 
same form as for the first-order coefficient ai(r; n). The only dependence on the order k 
in the BLM expansion is the constant accompanying the power of the renorm-group logs 



(inf^ + l)^ J (ln|^ + f ^ in Eqs. (6), (7) (details can be found in Ref. [5]). For the 
third order BLM, in particular, they were given in [7]: 
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Although the expressions for even higher orders in BLM is straightforward, in practice 
it is simpler to evaluate the BLM corrections in the widths directly in the running mass 
scheme. The general formalism is described in [5] , and using an additional technical trick 
given in Appendix simplifies necessary integrations. 

^Perturbative subseries including the coefficients of the form a„ = (^^^ 
are viewed as polynomials in n/. 
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2 Numerical estimates 

For the second-order non-BLM coefficient we have 



a2(r; a^) - af\r) = 5^ (r; fi) , Sr{0.25'; 1/4.6) -1.35 (9) 

non— BLM 

The simple interpolating expression accounting for the mass dependence of a^°^ , 



(10) 



5f (r;//) ~ -^zo(r) 1-4.5 (V^-0.25) + 13 (^^-0.25)^ + 0.7 -^--^ 
mf, L Vmft 4.D 

-12(V?-0.26)(A-i^)+71(v9-0.25f(A__L) 



is accurate enough in the relevant domain of quark masses and factorization scale /x.^ 

Using the evaluation 02°^ ~ 0.9 of the second-order non-BLM coefficients by Czarnecki 
and Melnikov [9] we obtain 

-non-BLM ^Q_252; 1 GcV) ~ -0.65 ^0(0.25^) (11) 

The second-order BLM coefficient has been addressed in the literature more than once; 
here the shift between /i — O and /i—lGeV constitutes about 0.78 /So for r~0.25^, and 

a^^^(r; /x) ~ -0.69 Pq Zo{r) at r = 0.25^ and £- = ^ (12) 

All BLM coefficients can be readily computed numerically along the lines of Ref. [5] 
using Eq. (A. 8) given below in Appendix. The required expression for the first-order 
perturbative correction with non- zero gluon mass is given in Ref. [10]. This yields 

af^(0.252;/.)~-0.52/32^o(r), af ^(0.25^; /.) ~ -0.27/5^ ^o(r) at (13) 

BLM-wise, it is not difficult to perform the all-order resummation of the perturbative 
corrections in terms of the running masses. The result reads 

^^^^(r;/.) ^ ^{^o(r) + ai(r;/.)^ + afLM(,.^^(a^)2 

+apM(^.^) (£^^%...| ^ 0.880; (14) 

it is independent of the scale for one starts with. 

The distribution of the gluon virtualities Q in rsi(6— >c) is illustrated in Fig. 1 in the 
case of pole masses (/x = 0) and for l GeV, /j, — 1.5 GeV and iJ, — 2 GeV. The effect of 
removing the infrared domain is self- manifest. Moreover, it is evident that a too significant 
part of the perturbative corrections with pole masses comes from gluon momenta below 
500 MeV, which would then cast legitimate doubts on the numerical results. 

^Recent analysis by DELPHI [8] yielded m6(lGeV)~4.6GeV and rwO.OS. 
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Figure 1: Gluon momentum scale distribution in T^ilb^c). Solid, dashed and dot-dashed lines 
correspond to n = lGeV, 1.5 GeV and 2GeV, respectively; lighter short-dashed line illustrates 
the case of /x = (pole masses) . The area under each curve gives the first-order perturbative 
coefficient. 



The growing component of the BLM series becomes sign-alternating starting the third 
or fourth order, this is a result of the infrared-free form of the width in the OPE-motivated 
approach. Consequently, the second-order approximation already turns out quite accu- 
rate, and the final number is practically reproduced by including the third-order term. It 
is also worth noting that using instead of af^{mb) a better physically justified expansion 
parameter [11, 3], for example the K-scheme a^^\mtj) or the 'dipole radiation' a^^\mij) 
yields a nearly precise perturbative factor already to the first order: 

A^^\r; At) ~ (1 - 0.105 - 0.018 + 0.017 - ...) . (15) 



The case of r^j'^'^ (6 — > u) is technically simpler. Here the full second-order corrections 
are known analytically [12], a2°'*(0) ~5.54, we have 52°'* — —6.9 and, therefore 

ai(0; 1 GeV) ~ -0.29 , a°°°-^^^(0; 1 GeV) ~ -1.35 , a^^^(0; 1 GeV) ~ 0.44 /3o ; (16) 

third- and fourth-order BLM corrections amount to 

af ^^(0; 1 GeV) ~ 1.23 , a^^^(0; 1 GeV) ~ 2.2 . (17) 

A meaningful BLM summation for b^u would require, however incorporating the four- 
fermion operator [13] 67^(1— 75)^ 'U7,^(l— 75)6 (5^i/— ^/iO^i^o) ^ well - it is strongly enhanced 
and has significant negative anomalous dimension. Accounting for this operator affects 
already the tree phase space coefficient at order ji^ /ml. 
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3 r„(B) 

In the complete OPE prediction for the semileptonic width, there remains a minor ar- 
bitrariness in incorporating the power corrections from the chromomagnetic and Darwin 
terms until their coefficients arc computed with the 0{as) accuracy.^ There are reasons 
to think that the factorized form has some advantages: 
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ml ml 



(18) 

where D{r) can be obtained from Ref. [14]: 

L>(r)=81nr + ^-^r-8r2 + ^r=^-^r^ ~ -18.3 (19) 

The kinetic and chromomagnetic expectation values here are finite-mQ matrix elements 
in actual B mesons. (As argued in Ref. [15], their shift compared to the limit mg— >oo is 
strongly suppressed.) 

To give the numerical prediction, we evaluate this at the canonical values mb(l GeV) = 
4.6 GeV, /^^(IGcV) = 0.4GeV^ /x^(lGcV) = 0.35 GcVl ^^^(IGeV) = 0.12 GeV^ and 
pi5(l GeV) = -0.15 GeV^ using r-0.252 and as{mb) = Q.22: 

l^c^l (tw)' (^)'(l-4.8[Br(B^X„£.)-0.0018]), ^ 0.0421 . 

(20) 

The above evaluation uses the full second-order perturbative corrections with BLM resum- 
mation Eq. (14), where the added non-BLM term is computed at q;s = 0.25. The resulting 

value of acquires extra factors 0.993 , 1.002 or 1.004 if one discards ©(a^) and 
higher terms, including only the third-order BLM correction, or both third and fourth 
orders, respectively. The dependence on the heavy quark parameters is as follows [16]: 

iKbl = X [1 + 0.4 (a, (m;,)- 0.22)] x 

[1 - 0.65 (m6(l GeV) -4.6 GeV) + 0.40 (mc(l GeV) -1.15 GeV) 
+ 0.01 (/x2-0.4GeV2) +0.10 (pf)-0.12GeV^) 

+ 0.05 (/x^-0.35 GeV^) - 0.01 (pi5+0.15 GeV^)] . (21) 



It should be noted that the Darwin expectation value p% appearing above is not a true 
matrix element normalized at 1 GeV, but is rather understood as the one extrapolated 
to jJL — Q: to two loops pf) p|,(l GeV) — 0.1 GeV^; the perturbative coefficients have 
been determined correspondingly. Modifying equations for using pf)(lGeV) proper is 

^The coefficient for the kinetic operator equals to the one in the parton term, while for the LS operator 
it is a combination of the latter and the coefficient for the chromomagnetic operator, see Ref. [2]. 
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straightforward. This does not change the value of \Vcb\ by any appreciable amount. The 
only drawback of the adopted simplified option is that the apparent extracted value of 
pI, is strongly correlated with the approximation applied. 

To summarize, the perturbative corrections to the total semileptonic widths of beauty 
particles are modest and allow precision theoretical control when viewed as the Wilson co- 
efficient of the leading heavy quark operator in the OPE. This assumes using well-defined 
low-scale running masses; these masses can be accurately extracted from experiment. As 
was pointed out long ago [3], the actual short-distance perturbative corrections to semilep- 
tonic widths are only about 10 percent, and can be accurately evaluated already at the 
one-loop level adopting a physically motivated perturbative coupling. 

Acknowledgments: 1 am grateful to M. Battaglia, D. Benson, I. Bigi, P. Gambino, 
T. Mannel and P. Roudeau for useful discussions. This work was supported in part by the 
NSF under grant number PHY-0087419. 



Appendix 

Here a brief description of technicalities involved is given following Ref . [5] . Computing 
the ELM corrections is most straightforward in the generahzed dispersive approach [17] 
employing the dispersion relation for the dressed gluon propagator: 



The form of the integrand suggests considering for a generic observable A the one-loop 
diagram with a non-zero gluon mass A 

A{X^)^l + ^A,{X') (A.2) 

TT 

and integrating it over A^ with the weight p(A^)^ = — Im q;s(— A^) 4^: 
A simple trick appears useful. Expressing 

with an arbitrary M, the integral of the first term yields identically Ai{0)^^^^-^ and, 
therefore 

= l + MO)^-l f p(A^) (a,(0)-^-A,(A^)) . (A.4) 
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The standard BLM expansion amounts to using literal one- loop q;s(/c^): 



l + *^ln^ ■ 



(A.5) 



As had been shown in [18], in this setting the dispersive approach reproduces the literal 
BLM series. 

Taking M'^ — es niQ we obtain expansion directly in terms of af^{mQ) : 
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where as = af^{mQ), and 



(A.6) 



(A.7) 



This form is convenient for numerical integration. It also allows the straightforward 
expansion in asimq): 



iBLM 



n=0 

n 

2fe+i fdX^ 



.n+2 



ik=0 •' ^ 
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A2 3 



n-2k 



MO)^^^^.-Am]- (A.8) 



Since the infrared part is removed from the one-loop diagram for the total width, the 
corresponding Ai{\'^) is a real analytic function in the vicinity of zero, and no ambiguity 
in Eq. (A.6) is encountered. Nevertheless, as follows from the generating integral (A.6) 
the series (A.8) is only asymptotic. 

Computing BLM corrections with the running masses toq(/x) requires expressions for 
the order-Q;^ terms in A ^""^ (^) and fi'^'"""' (fi) . They are given by integrating the SV spectral 
density (Eqs. (16), (19) of Ref. [5]) with the proper power of energy: 



A-'(/.; A) = 1^^(;,2-A^) 



1 — ^) y iJi? — X^ — ^ A + 7 A arctan 



A 



8 



A) 



4ct, 
37r 



(A.9) 



As discussed in Refs. [18, 17], the quantity _;^2 dA]^(A ) -j^^ qualitatively viewed 
as describing the contribution of gluons with virtuality A^ in the one-loop process. This 
convention has been adopted in Fig. 1. 
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